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. . . Abstract. We review computational approaches to understanding the origin of the Initial Mass 

, Function (IMF) during the formation of star clusters. We examine the role of turbulence, gravity 

• and accretion, equations of state, and magnetic fields in producing the distribution of core masses 

' - the Core Mass Function (CMF). Observations show that the CMF is similar in form to the 

, IMF. We focus on feedback processes such as stellar dynamics, radiation, and outflows can 

O ' reduce the accreted mass to give rise to the IMF. Numerical work suggests that filamentary 

^ accretion may play a key role in the origin of the IMF. 
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1. Introduction 

The Initial Mass Function (IMF) plays a central role in astrophysics because it encap- 
sulates the complex physics of star formation. Observations suggest that the IMF can be 
vario usly described as a piece- wise power- law ( Kroupall2002l ). a logn ormal ( Miller &: Scalo 
or a lognormal distributio n with a high-m ass power- law tail (|Chabrieill2003l) . The 
high mass behaviour of the IMF ISalpeteii (|l953f) is a power-law, dN cx m~'^'^dm for stel- 
lar masses m ^ 0. 5Mndnt- The pe ak mass for isolated stars in the galactic disk is 0.1 Af0 
and 0.2 - 0.3Mq (|Chabrieij|2003t ) for the bulge. The form of the IMF is similar in many 
different galactic and extragalactic environments such as globular clusters, where in one 
has a large range of metallicities and concentrations ( Paresce fc De Marchi l2000l ). The 
current evidence therefore tends to support the notion that the IMF is universal. 

What physical processes produce the IMF? Obser vations show that it emerges dur- 
ing th e early stages of the formation of star clusters ( Meyer et al. 2000t Zinnecker et al 



19931) . Young stars are formed within gravitationally bound subunits of a cluster- forming 



environment known as "cores" whose mass distribution - the core mass function (CMF) 
- strongly resembles the IMF. Cores are closely associated with fila ments, as recent ob - 
servations using the Spitzer and Herschel observatories clearly show ( Andre et al.ll2010l ). 
Competing physical processes such as turbulence, gravity, cooling and thermodynamics, 
as well as magnetic fields play significant roles in building the CMF, as well as filaments, 
within molecular clouds. Feedback processes such as radiation from massive stars, jets 
and outfiows, as well as stellar dynamics serve to truncate the accretion of material onto 
stars and their natal disks. These may lead to the emergence of the IMF from the CMF. 

This review focuses on the critical role that computation is playing in exploring the 
origin of the IMF in clusters. We focus first on the processes leading to the CMF, and 
then discuss feedback processes that may convert the CMF i nto the IMF. Recent r eview s 
of the computationa l aspects of th e IMF may be found in Mac Low fc KlessenI (j2004l ). 
Bonnell et all (|2007l ). lLarsonl (|2007t ). and lKlessen et all (l2009l ) The theory of the IMF is 



covered by Hennebelle (this volume), and McKee fc Ostrikei ( 2007t ) 
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2. Structure formation and the IMF 

Structure formation in the diffuse ISM as weh as the dense molecular medium is a 
shocking affair. Diffuse atomic hydrogen near the midplane of our Milky Way is ob- 
served to be organi zed as a plethora o f filamentary structure, bubbles, supernova rem- 
nants, HII regions (jTavlor et al. I I2OO3I) . Filamentary structure also characterizes giant 
molecular clo uds - as seen in the large scale extinction maps of the Orion and Mono- 
cerus clouds ( Cambresv 1999() . On smaller scales, a 850 micron continuum map of a 
lOpc region in Orion sho ws that Bonner-Ebert like cores are associated with filaments 
(| Johnstone fc Ballvl 120061 ). Herschel observations, such as those of the IRDC filament 



known as the " snake" , show massive stars and a star cluster in formation (jHenning et al 
I2OIO). 

How can cloud structure be characterized? One fruitful approach is to measure the 
probability distribution functions (PDFs) of the cloud column density of all of the gas 
in a given molecular cloud which is readily measured directly from the extinction data. 
This PDF for clouds without star formation (such as Lupus V and the Coal Sack clouds) 
turns out to be lognormal, wh ereas that of sta r forming clouds is a lognormal plus high 
mass, power-law tail (jKainulainen et al. I I2OO9I) . This is interpreted as arising from the 
effects of gravity which drives collapse. 

Filaments have various origins and arise in supersonic turbulent media due to the 
intersection of shocks waves, by gravitational break-up of self-gravitating sheets, or by 
thermal instabilities of various kinds. A physically plausible picture for the origin of the 
CMF is that filaments with sufficiently high values of their mass per unit length produce 
cores by grav itational instabifity. Gravitational fragm entation of filaments is well studied 
theoretically ([Nagasaw al llQSTt iFiege fc Pudritdlioom but has received renewed emphasis 
with the Herschel observations. 

The similarity in the functional form of the CMF and th e IMF has been observed in 
many clouds starting with the study of iMotte et al.l (|l998[ ) in p Oph. The CMF in the 



Pipe Nebula, as another ex ample, can be shif ted into the IMF by converting a fraction 
e ~ 1/3 of its mass to stars (jAlves et al.ll2007l ). 



3. From clouds to the CMF 



3.1. Turbulence 

Supersonic turbulence rapidly compresses gas into a hierarchy of sh eets and filaments 
wherein the denser gas undergoes gravitational collapse to form stars (IPorter et al. 1994: 



Vazquez-Semadeni et al.ll995t Ostriker et al.ll999t Klessen fc Burkerti200ll : Padoan et al 



2OO1I iBonnell et al.ll2003l : iTiUev fc Pudritzll2Q04t iKrumholz et al.ll2007l ). Supersonic gas 
is also highly dissipative and without constant replenishment, damps within a crossing 
time. There are many sources of turbulent motions that can affect molecular clouds, 
such as galactic spiral shocks in which most giant molecular clouds form, supernovae, 
expanding HII regions, radiation pressure, cosmic ray streaming, Kelvin-Helmholtz and 
Ray leigh- Taylor i nstabilities, gravitational in stabilities, and bipolar outflows from regions 
of star formation ( Elmegreen fc Scald 2004 ). 

Supernova driven bubbles and turbulence in the galac t ic dis k have been simulated by 
several groups. The work of Ide Avillez fc Breitschwerdt (2004) simulates the high reso- 
lution (down to 1.5 pc scales) global structure of the ISM as a function of the supernova 
rate. Densities range over six orders of magnitude, 10"** ^ n ^ 10^ cm"'^ and multi- 
phase density PDFs are given. Simulations of the global structure of a supernova lashed. 
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multi-component, instellar medium ([Tasker &: Brvanll2008l : IWada fc NormanI 120071 ) find 
that the density PDF follows a lognormal distribution. 

Supersonic turbul ence produces hierarchica l structure that can be described by a log- 
normal distribution (|Vazquez- Semadeni|[l99l . A lognormal arises whenever the proba- 
bility density of each new step in density increment in the turbulence is independent of 
the previous one. As an example, consider a medi um that undergoes a series of random 
shocks whose strengths are uniformly distributed ( Kevlahan fc Pudrit j 20091) . The den- 
sity at any point is the product of the shock-induced, density jumps. Taking the log of 
this relation, and then the limit of a large number of shocks, the central limit theorem 
then shows that the log of the gas density should be normally distributed - hence the 
lognormal. It can also be shown mathematically that the convergence of the distribution 
for a finite number of shocks is very rapid - just 3 or 4 shocks will give a distribution 
that is very highly converged to a lognormal. 

Lognormal behaviour for the CMF has been found in a wide variety of simulations using 
various types of codes (eg. SPH, AMR) and setups (driven or not driven, periodic boxes, 
initial uniform spheres, etc.). Early results showed that lognormal behavi our in periodi c 
box simulations is independent of details on how the turbulence is driven (lKlessenll200l . 

It is interesting that lognormal distributions appear across science (eg. physics, bi- 
ology, medicine, etc.), not just in fluid mechanics. The key difference between normal 
and lognormal distributions in general is that the former arises for additive processes, 
whereas the latter arise in multiplicative ones (jLimpert et al ] l200ll) . As a concrete ex- 
ample, consider a simple dice game where one first adds the values on the faces of two 
thrown dice - the distribution of results (ranging from 2 to 12) is a normal distribution 
whose mean is 7. If one multiplies the two values however, the resulting distribution of 
numbers is highly skewed (ranging between 1 and 36) and is described by a lognormal. 

Thus, lognormals characterize structure in the diffuse ISM as well as in denser molecu- 
lar gas because shocks are the dominant process for configuring the gas. This is indepen- 
dent of exactly how the shocks are produced. Of physical significance are the mass of the 
peak of this distribution, and its width (Tq. The latter depends on both the thermal state 
of the gas as well as the rms Mach number of the turbulence. The standard deviation ctq 
is found from the simulations, and takes the form; = ln{l + b'^AP ) where M is the rrn s 
Mach number of the turbulence and b ~ 0.5 is a fitting parameter (|Padoan et al.lll997 ). 



3.2. Gravity and accretion 

Adding gravity to turbulence changes th is distribution - a power-law tail app ears at 
the high mass end of the si mulated CMF (iLi et al.ll2003t iTillev fc Pudrit jl2004| ). In the 
semi-analytic treatment of Padoan fc Nordlund J 2002fl . the power- law a rises because of 
the turbulence spectrum of the turbulence. In Hennebelle fc Chabrier ( 20081) . a Press- 
Schecter formalism is adopted to argue that a lognormal plus power-law behaviour is the 
consequence of imposing a star formation threshold (eg. Jeans' criterion) and gravity at 
high mass. 

Much of the debate concerning the origin of the CMF a nd IMF has focused on whether 
stars form by co mpetitive accretio n ( Bonnell et al.ll200ir) or by the gravitational collapse 
of discrete cores ( Krumholz et al. 2005f) . Cluster formation simulations often use initial 
top-hat density profiles (uniform spheres) that are chosen to mimic the observed initial 
conditions of cluster forming clumps (eg. a hundred solar masses of material, at temper- 
ature of lOK, size of half a pc, and mean density of at least 10^ cm^"^). Such simulated 
clumps start t o undergo global gravita.tiona l collapse in less than a free-fall time (eg. 
Bonnell et al] (|200l[ ). [Tiilev fc Pudrit3 |2004l )) as the turbulent energy is dissipated. The 
collapsing background ramps up the density of the gas including those in fiuctuations 
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such as the filaments. This drives up the accretion of gas into the filaments pushing some 
of them towards gravitational fragmentation. Sink-particles (taken as proxies for cores) 
first appear within the filaments. Cores will collapse more quickly than the collapsing 
background clump because they are denser and therefore have shorter free-fall times. The 
the collapse carries the collection of sink particles and filaments into the ever deepening 
potential well of the clump. 

Cores are not isolated objects within filaments. Rather, they continue to undergo 
considerable filamentary accretion. The first o bjects to appea r with in most simulations 
generally become the most massive cores. In Baneriee et al.l ( 20061 ). the rapid growth 
of the first, and most massive star by filamentary accretion was followed with an AMR 
(FLASH) code. These simulations showed that the filament was formed at the intersection 
of two sheets- shocks. 



3.3. Equations of state 

The thermal state of the bulk gas is a major factor in determining the local Jeans mass 
from point to point in molecular clouds. Both theory and simulations show that equation 
of state plays a very important role in controling the gravitational fragmentation of 
the gas. The local Jeans mass scales with the local temperature and density as AIj oc 
7^3/2^1/2 ppj. sijjjpig polytropic equations of state P (x , the Jean's mass can be written 
purely as a power law of the density; Mj oc p3/2(7~(4/3))^ This scaling suggests that for 
7 > 4/3, the Jeans mass increases with density, which puts an end to fragmentation. 
Indeed, simulations show that strong fragmentatio n prevails for 7 ^ 1, gets progressively 
weaker for 7 > 1, and stops altogether for 7 > 1.4 (jLi et al ] |2003l) . 

Local energy sources can raise the fragmentation mass by changing the temperature of 
the region. Thus, the accretion luminosity released by massive stars in particular, must 
reduce the degree of fragmentation within a l ocalized region a,round such a heat source, 
as has been demonstrated by several groups ([Krumholz et al. I l2007h . For massive stars, 
this region is limited in extent - roughly 1000 AU or so (see §4). 



3.4. Magnetic fields 

Magnetic fields play several different roles in star formation. As has long been known, 
they can control the gravitational stability of a gas if their energy density exceeds that 
of gravity. This is formalized by mass to flux ratio; F = 2'k\/GYj/ B — lA/3^/^nJ'^ where 
nj is the number of Jeans masses and (3 is the ratio of gas to magnetic pressure. The 
fragmentation of a uniform cloud is highly suppressed for subcritical clouds (F < 1). 
Simulations of slightly supercritical clumps in uniform magnetic fields show that the 
field channels the collapse into large sheets that are perpendicular to the direction of 
the field. Slightly m ore supercritical clouds however, break up into more substructure 
including filaments ( Tillev fc Pudritzl 2007 ) . Turbulence creates a very broad range of 
magnetizations of cores. This is because shocks sweep material along field lines where it 
accumulates in filaments - increasing the mass to flux ratio in those regions, and greatly 
reducing the mass to flux ratio in the more diffuse zones left behind. Simulations show 
that initially supercrit ical magnetized clouds results in cores that range from critical to 
strongly supercritical (Padoan fc N ordlund 2002 : Tillev fc Pudritz 2007 ). in agreement 
with the observations ( Crutchcr 2 0071 ). This may also account for the fact that magnetic 
fields are more dominant in the diffuse gas than in molecular gas. 

When turbulence is added to a subcricital cloud, the column density PDF of a cloud is a 
lognormal but with a very small stand ard deviation. When am bipolar diffusion is added, 
the lognormal broadens considerably ( Nakamura fc Li 20081 )). Supercritical subregions 
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can form within the subcritical cloud, and it is within these regions that star formation 
can proceed, albeit at a very heavily reduced rate. 

The second major aspect of magnetic fields is that twisted fields exert torques on gas 
and can therefore transport angular momentum away from spinning bodies. The origin 
of an gular momentum in supersonic turbulence is tha t oblique shocks produce spinning 
cores ( Jappsen fc Klessenll2004 : Tillev fc Pudrit dr2007 )). Turbulence simulations produce 
a broad distribution of angular momenta of cores ranging over nearly two orders of 
magnitude. This implies that a broad distribution of disk sizes should result from the 
collapse of these systems - from very small to very large disks. The distribution of angular 
momentum vectors is also quite varied and is not particularly aligned with the filament 
principal axes. 

Magnetic fields participate in several kinds of braking. On the level of the cores, tor- 
sional Alfven waves h ave long been known to be able to extract significant amounts 
of angular momentum Basu fc Mouschoviad (1993). The collapse of rotating cores pro- 
duces ina£netized_outflows as demonstrated in a variety of initial co re models; cylin- 
ders (lTomisakall2002l). Bonner -Ebert spheres (|Baneriee fc Pudritj|2006l). uniform s pheres 
( Hennebelle fc Fromang|2008 ) , and singular isothermal spheres ( Mellon fc Lil2008l) . Early 
outflows could sweep up significant amounts of material and have been implicated as the 
basic physics in the CMF to IMF efficiency factor e (see §4.3). 

Finally, the combination of radiative and MHD effects (RMHD) limits the fragmenta- 
tio n of cluster-for i ning g as. Attempts to model RMHD on cluster scales have been made 
bv IPrice fc Bate ( 2Q09t ) who used Euler potentials to approximate the MHD in SPH. 
This work shows that the MHD in their code strongly supresses fragmentation even for 
clouds that are fairly supercritical (eg. F = 3). On smaller scales, grid-based RMHD 
methods applied to collapse and outflows show that radiative heati ng (in fiux limite d 



diffusion limit) makes su bstantial changes to the extent of outflows (jCommeryon et al 
2010HTomida et al.ll20ld ). 



4. Feedback: from the CMF to the IMF 



4.1. Stellar dynamics and filamentary accretion 

The collapse of lO^Af© clumps pulls both the material in the filaments and the sink 
particles deeper into the central potential well. The sink-sink interactions, which are 
modeled to be N-body gravity then come into play and are responsible for the creation 
of a small stellar cluster. The competitive accretion scenario predicts that these objects 
comp ete for gas in the de nse centre of the gravitational potential well, resulting in an 
IMF (jBonnell et al.l 120011 ). For larger clumps (10^ — IO'^Mq) the evolution of a cluster 
form ing region is a h ierarchical process in which subclusters form and ultimately merge 
(B onnell et al.ll2003l ). The resulting stellar interactions are frequent and close enough to 
truncate protostellar disks. 

Sink particles as implemented by IFederrath et al. I (l2010l) have the additional desirable 
feature that they form in local potential minima giving them a more hydrodynamic 
character. This has an important consequence. As opposed to competitive accretion onto 
particles in a general potential well, filamentary accretion largely ends when the N-body 
interactions become strong enough to kick the sink particles ou t of their feed i ng filam ents 
Evidence for filamentary accretion is seen in the simulations of Duffin et al.l ()2010al ) , and 
shown in Figure 1 (featuring the collapse of a lOOM© ~ lOOMj initial tophat clump). The 
first frame, shows that sink particles form in filaments. The second shows that dynamical 
interactions have started between them as the collapse of the clump proceeds, and the 
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Figure 1. Cluster formation and filamentary accretion. Top left: sink particles form in filaments. 
Top right: sink particles begin dynamical interaction. Bottom left: cluster dynamics, accretion 
has ceased. Bottom right: ma ss accretion history for each sink particle - accretion ends when 
dynamical interactions begin (|Duflin et al.ll2010al ) 



third shows the dynamical end state. The final frame shows that the accretion histories 
of each particle shuts off when when dynamical interactions become important. 

4.2. Radiative feedback 

Comprehensive 2D simulations of massive star f ormation including fu ll frequency, radia- 
tive feedback effects and dust are presented by lYorke fc Sonnhalteil ([2002). For cluster 
formation, a 3 D treatment of radiati ve transfer in a highly inhomogenous medium is 
essential. Since iKrumholz et al.l (|2007l ). it has become widely appreciated that radiative 
heating of the gas is needed to prevent excessive gravitational fragmentation. Radiative 
heating of the gas prevents the filaments from fragmenting as much as they might. Their 
gas drains primarily into the central massive forming disk and star, and fragmentation 
out to 1000 AU scales is prevented. These calculations still invoke grey, flux limited dif- 
fusion so there is still a need to ex amine the role of different frequency regimes in this 
process. Simulations by iBatd (l2009l) confirm this picture using SPH techniques - wherein 
suppression of brown dwarfs by a factor of 4 is often observed compared to simulations 
without radiative feedback. The ionizing radiation from massive young stars becomes 
important during cluster formation as HII regions drive hot ioinized fiows in the cluster 
environment. Results show that companions t o massive stars limit the growth of the 
latter by a starvation effect ( Peters et al. 2010l) . 
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4.3. Outflows 

Outflows carry substantial amounts of mechanical energy which, if coupled into clump dy- 
namics ^_couMafiectJheJMF. Two aspects of feedback are important outflows in a cluster 
region ("Kless en et al. 20091 ). The first is at the core level, in which an outflow removes the 
collapsing gas at some efficiency. If this efficiency is quite high (so t hat e :^ 1/3), then 
one m ay resolve the hypothesized conversion of the CMF to IM F dMatzner fc McKee 
20001) . Simulations show that the efficiency may be much lower (jPuffin et al.l bOlOb ). 



The second is at the level of the clump and the question as to whether or not the collec- 
tion of outflows can co ntinue to excite turbulence and thereby regulate cluster formation 
Norman fc silS (|l980l) . 



The situation at the local level is still unclear. Simulations bv iBaneriee et al.l (|2007f ) 
show that super s onic t urbulence is not driven by magnetized jets. At the level of clumps, 
Nakamura fc Lil ()2008l ) included ambipolar diffusion and outflows in a cloud that is nearly 



critical (F = 1.1 initially), and found, not unexpectedly, t hat the magnetic fie lds regulated 
a rather slow rate of star formation. On the other hand, Wang et al. ( 2010l ) used a mass 
to flux of F = 1.4 in a code without AD, and observed long time regulation of the cluster 
formation by outflows which maintained turbulence. These results may all depend upon 
the limited modeling of the the full dynamics of outflows that have been incorporated 
into the simulations. 



5. Synthesis: the IMF of clusters 

Numerical simulations have become the primary tool with which to investigate the 
origin of the IMF in clusters. Turbulence, the equation of state, and gravity play the key 
roles of fllamenting the gas and breaking it into accreting cores describable by a lognormal 
distribution with high mass power-law tail. Radiative feedback controls fragmentation 
rates at the low mass end of the IMF whereas feedback by outflows may not be as efficient 
as previously claimed in converting the CMF to the IMF. Filamentary structure from 
clouds to clumps may turn out to play a key role in the entire process. 
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